Charged particle's elliptic flow in 2+1D viscous hydrodynamics at LHC ( \fs= 2.76 

TeV) energy in Pb+Pb collision. 
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In the Israel-Stewart's theory of second order hydrodynamics, we have simulated yfs=2. 76 TeV 
Pb+Pb collisions. ALICE data for the centrality dependence of charged particles multiplicity, 
Pt spectra in 0-5% collisions, centrality dependence of integrated and differential elliptic flow are 
analysed. Analysis indicate that while ALICE data on charged particles multiplicity or pr spectra 
do not demand any viscosity, viscosity is demanded by the elliptic flow data. From a simulataneous 
fit to all the data sets, viscosity to entropy ratio in ^/s NN =2.76 TeV Pb+Pb collisions is extracted 
as, n/s = 0.06 ± 0.02. 

PACS numbers: 12.38.Mh ,47.75.+f, 25.75.Ld 



I. INTRODUCTION 

Elliptic flow is one of the most important observables 
in relativistic energy heavy ion collisions. Elliptic flow 
measure the azimuthal correlation of produced particles 
with respect to the reaction plane. Finite elliptic flow in a 
non-central collision is now regarded as a definitive signa- 
ture of collective effect in relativistic energy nuclear col- 
lisions P, Q . Qualitatively, elliptic flow is best explained 
in a relativistic hydrodynamical model, rescattering of 
secondaries generates pressure and drives the subsequent 
collective motion. In a non-central collision, the reac- 
tion zone is asymmetric (almond shaped), and pressure 
gradient is large in one direction and small in the other. 
The asymmetric pressure gradient generates the elliptic 
flow. As the fluid evolve and expands, asymmetry in the 
reaction zone decreases and a stage arise when reaction 
zone become symmetric and system no longer generate 
elliptic flow. Elliptic flow is an early time phenomena. It 
is a sensitive probe to, (i) degree of thermalisation, (ii) 
transport coefficients and (iii) equation of state of the 
early stage of the fluid 0, H, 0] ■ Relativistic hydrody- 
namics has been successfully applied to analyse various 
experimental data in Au+Au collisions at RHIC energy 
(i/s=200 GeV). Hydrodynamical evolution of near ideal 
QGP fluid, with shear viscosity to entropy ratio rj/s wO.l, 
initialized to central energy density Si f=s 30 GeV/ fm 3 at 
initial time Tj «0.6 fm/c, is consistent with experimen- 
tal elliptic flow in central and mid-central collisions [H-Q . 
Small viscosity is consistent with string theoretical cal- 
culation based on ADS/CFT correspondence [l(| that 
viscosity of a strongly interacting matter is bounded by 
the KSS bound rj/s > 1/4tt. One wonders about the 
viscosity to entropy ratio {rj/s) at LHC energy. Present 
paradigm is that rj/s has a minimum, possibly with a 
cusp, around the critical temperature T = T c Since, 
at LHC, produced fluid will be at a higher temperature 



* E-mail: victor@veccal.emet. in 
tE-mail akc@veccal.ernet.in 



than the fluid at RHIC, one expects higher value of rj/s 
at LHC. 

Recently, ALICE collaboration measured charged par- 
ticles elliptic flow in 10-50% Pb+Pb collisions at LHC 
energy, v /s Arjv =2.76 TeV [llj]. Differential elliptic flow in 
10-50% collisions at LHC is also rather similar to that 
at RHIC. Integrated elliptic flow on the other hand in- 
creased by ^30%. ALICE collaboration also measured 
charged particles multiplicity [l]|, pr spectra [l4| in 
Pb+Pb collisions. At LHC multiplicity increase by a fac- 
tor of ~2.2, though the centrality dependence is rather 
similar to that in -^^^=200 GeV Au+Au collisions. 
High px suppression is also more at LHC energy. In 
0-5% collisions, pt spectra of charged particle's are sup- 
pressed by a factor Raa ~ 0.14 at p T —6-7 GeV, which 
is smaller than at lower energies. In peripheral collision, 
suppression is modest, Raa ~0.6-0.7. 

Increase of elliptic flow at LHC energy collisions has 
been predicted in several hydrodynamical and hybrid 
models. In [l5l - fT7j | it was predicted that flow will increase 
by 20-30%. Some models also predicted decreased ellip- 
tic flow in LHC energy [I8I - I20I ] . For example, in [2(| , one 
of the present author predicted that at y/s NN =5 TeV 
Pb+Pb collisions, elliptic flow will decrease by ~15%. 
The prediction was based on several assumptions, which 
are incorrect as of now. For example, contrary to the 
present understanding that multiplicity increases with 
energy by a power law, in (20| a logarithmic dependence 
was assumed. This led to much reduced multiplicity 
than observed experimentally. Also, in [20| . confinement- 
deconfinement transition was assumed to be 1st order, 
contrary to the present understanding that the transi- 
tion is a cross-over. 

In the present paper, in a viscous hydrodynamic model, 
we have analysed the ALICE data for the centrality de- 
pendence of the charged particles multiplicity, px spectra 
in 0-5% collisions, centrality dependence of integrated 
and differential elliptic flow and extracted viscosity to 
entropy ratio for the fluid produced in ^/s NN =2.1Q TeV 
Pb+Pb collisions. Analysis indicate that the ALICE data 
on the centrality dependence of charged particles multi- 
plicity and 0-5% px spectra are best explained in ideal 
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hydrodynamics. Elliptic flow data on the otherhand de- 
mand viscosity. Small viscosity rj/s = 0.06 ± 0.02 is re- 
quired if all the data sets are analysed simultaneously. 

The paper is organised as follows: in section [Til we 
briefly explain the second order Israel-Stewarts theory 
of dissipative hydrodynamics. Hydrodynamics equations 
are closed with an equation of state (EoS). In section ITO1 
we discuss the EoS used in the present study. Section ITVl 
deals with the initial conditions used for hydrodynamic 
simulation of v / s JVA ,=2.76 TeV Pb+Pb collisions. Simu- 
lation results are compared with experiments in section 
IVl Finally in section fVTl we summarise the results. 



II. HYDRODYNAMIC MODEL 

In Israel-Stewart's second order theory (2lj . thermody- 
namic space is extended to include the dissipative fluxes. 
Relaxation equations for the dissipative fluxes are ob- 
tained from the entropy law, <9 M s M > 0. Equation of 
motion of the fluid is then obtained by solving simultane- 
ously the conservation equations, e.g. energy-momentum 
and the relaxation equations for the dissipative fluxes. In 
^^^=2.76 TeV Pb+Pb collisions, the central rapidity 
region is approximately baryon free and dissipative ef- 
fect like conductivity can be neglected. We also neglect 
the bulk viscosity. In general bulk viscosity is much less 
than the shear viscosity. But in QCD, near the transi- 
tion point bulk viscosity can be large |22|. [23| . Effect of 
bulk viscosity on elliptic was studied in [2J] . The effect is 
not large. We have also extended the code AZHYDRO- 
KOLKATA to include bulk viscosity. Preliminary re- 
sults are reported in (2||. We also find that inclusion 
of bulk viscosity marginally modify the spectra or ellip- 
tic flow. Equation of motion of the fluid in y^s NN =2. 76 
TeV Pb+Pb collisions is then obtained by solving, 



Dtt^ 



= 0, 



-{it> tv - 2rjV 



[u^ vX + u v ^ X ]Du x . 



(1) 



(2) 



EqfT] is the conservation equation for the energy- 
momentum tensor, = (e + p)u^u v — pg^ v + n^' y , 
e, p and u being the energy density, pressure and fluid 
velocity respectively. ■k^ v is the shear stress tensor (we 
have neglected bulk viscosity and heat conduction) . Eq|2] 
is the relaxation equation for the shear stress tensor . 
In EqfJl D = u^d^ is the convective time derivative, 
y<» u v> = i + y^) - ±(d.u)(g^ - u^u") is a 
symmetric traceless tensor. r\ is the shear viscosity and 
T n is the relaxation time. It may be mentioned that in 
a conformally symmetric fluid relaxation equation can 
contain additional terms [24| . 

Assuming boost-invariance, EqsfT] and [5] are solved in 
(t = \Jt 2 — z 2 , x, y, rj s = ^lnf+zl) coordinates, with the 
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FIG. 1: (color online) Black and red circles are Wuppertal- 
Budapest lattice simulations [2£| for entropy density and pres- 
sure. Solid lines are the fit to the lattice simulation with 
model the EOS (Eq[3}. For completeness energy density in 
the model EOS is also shown. 



code '"AZHYDRO-KOLKATA", developed at the Cy- 
clotron Centre, Kolkata. Details of the code can be found 
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III. EQUATION OF STATE 

Hydrodynamic equations are closed with an equation 
of state p = p(e). Most reliable information about the 
QCD equation of state is obtained from lattice sim- 
ulations. It is now established that the confinement- 
deconfinement transition is not a thermodynamic phase 
transition, rather a cross-over (26[. Since it is cross-over, 
there is no unambiguous temperature where the transi- 
tion take place. Inflection point of the Polyakov loop 
is generally quoted as the (pseudo) critical temperature 
for the confinement-deconfinement transition. Currently, 
there is consensus that pseudo critical temperature for 
confinement-deconfinement phase transition is T c «170 
MeV H3-[l|. 

We have constructed EoS by complementing 
Wuppertal-Budapest lattice simulations [28| with a 
hadron resonances gas comprising all the resonances 
upto mass m res « 2.5 GeV. Entropy density is parame- 
terised as, 



s = 0.5[(1 - ta,-nh(x)]s H RG + 0.5[1 + tanh(a;)]s/ attJce (3) 

where srrg and Siattice are entropy density of the con- 
fined and deconfined phase, x — , AT C = 0.1T C . 
From the parametric form of the entropy density, pres- 
sure and energy density can be obtained using the ther- 
modynamic relations, 



P(T) 



s(T')dT' 



e{T) = Ts-p. 



(4) 
(5) 
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As stated earlier for srrg we have used hadronic 
resonance gas comprising all the resonances with mass 
m res < 2.5 GeV. For siattice we use a parametric form 
for the Wuppertal-Budapest simulations for entropy den- 
sity, 



S lattice 



[a + f3T] [1 + tanh 



T-T c 
8 



(6) 



with cv=6.74, 0=0.32, T c =174 MeV and 6=70 MeV. 

In FigHJ Wuppertal-Budapest simulations for the en- 
tropy density and pressure, as a function of temperature 
are shown. The solid lines are fit obtained to the lat- 
tice simulations for entropy density and pressure in the 
model EoS. The model EoS correctly reproduces lattice 
simulations for entropy and pressure. 



IV. INITIALISATION OF THE FLUID 

Solution of partial differential equations (Eqs lll2j) re- 
quires initial conditions, e.g. transverse profile of the 
energy density (s(x,y)), fluid velocity (v x (x,y),v y (x,y)) 
and shear stress tensor (n^ (x,y)) at the initial time Tj. 
One also need to specify the viscosity (rj) and the relax- 
ation time (t„). A freeze-out prescription is also needed 
to convert the information about fluid energy density 
and velocity to particle spectra and compare with ex- 
periments. 

In a recent publication, we have analysed the cen- 
trality dependence of charged particles multiplicity in 
y/s NN =2.76 TeV Pb+Pb collisions [30(. It was shown 
that very short thermalisation time Ti=0.2 fm or very 
large thermalisation time r^=4 fm, is not consistent with 
the ALICE measurements. Presently, we assume that the 
fluid is thermalised at the time scale Ti=0.6 fm. Initial 
fluid velocity is assumed to be zero, v x (x, y) = v y (x, y) = 
0. At the initial time Tj, in an impact parameter b col- 
lisions, we distribute the initial energy density as in a 
Glauber model 



e(b,x,y) = e l [{l-x)N part (h,x,y)+xN coU (b,x,y)}, (7) 



£j in EqJ7]is the central energy density in impact param- 
eter b = collision. N part and N co u are the transverse 
profile of the average number of participants and aver- 
age number binary collisions respectively, calculated in 
a Glauber model, x is the fraction of hard scattering 
and is an important parameter for elliptic flow develop- 
ment. Initial spatial eccentricity is more if initial density 
is dominated by the hard scattering, so is the elliptic 
flow [31j . In Glauber model of initialisation with energy 
density, unless the hard scattering fraction is large, AL- 
ICE data on the centrality dependence of the charged 
particles multiplicity is not reproduced. We assume 90% 
hard scattering fraction. The shear stress tensor was ini- 
tialised to boost-invariant values, ir xx = tt vv = 2r]/3Ti, 
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FIG. 2: ALICE measurements for centrality dependence of 
charged particles multiplicity are compared with hydrody- 
namical simulations for t;/s=0-0.16. 



tt xv =0. For the relaxation time, we used the Boltzmann 
estimate tv = 3ij/2p. Finally, the freeze-out temperature 
was fixed at Tp=130 MeV. Elliptic flow is an early time 
phenomena and depend little on the freeze-out tempera- 
ture. In viscous hydrodynamics, particle production has 
contributions from the non-equilibrium part of the dis- 
tribution function. We have checked that at the freeze- 
out, non-equilibrium contribution to particle production 
is much smaller than the equilibrium contribution. We 
note that the freeze-out condition does not account for 
the validity condition for viscous hydrodynamics, i.e. re- 
laxation time for dissipative fluxes are much greater than 
the inverse of the expansion rate, trO^u^ << 1. Re- 
cently, Dusling and Teaney 32j implemented dynami- 
cal freeze-out condition. In dynamical freeze-out, non- 
equilibrium effects are stronger than obtained at fixed 
temperature freeze-out. 



TABLE I: Central energy density (e,) and temperature (Ti) at 
the initial time n=0.6 fm/c, for different values of viscosity 
to entropy ratio (rj/s). The bracketed values are estimated 
central energy density and temperature in y / s JVJV =200 GeV 
Au+Au collisions [9j. 








0.08 


0.12 


0.16 


Si 


143.0 ± 6.0 


126.0 ±5.5 


115.0 ±5.2 


103.0 ±4.6 




(35.5 ± 5.0) 


(29.1 ±3.6) 


(25.6 ±4.0) 


(20.8 ± 2.7) 




548 ±5 


531 ±7 


520 ±6 


504 ± 7 


(MeV) 


(377 ± 14) 


(359 ± 12) 


(348 ± 14) 


(331 ±11) 



The only parameters remain to be fixed is the central 
energy density and viscosity coefficient rj. We assume 
that the ratio of viscosity to entropy tj/s remain a con- 
stant through out the evolution and simulate Pb±Pb col- 
lisions for four values of viscosity, (i) r//s=Q (ideal fluid) 
(ii) rj/s = 1/4tt « 0.08 (KSS bound), (iii)r?/s=0.12 and 
(iv) ?y/s=0.16. During the evolution, fluid cross over 
from QGP phase to hadronic phase. The assumption 
that rj/s remain unchanged during the evolution vio- 
lates the present understanding that rj/s has a minimum 
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FIG. 3: ALICE measurements for charged particles pr dis- 
tribution 0-5% Pb+Pb collisions are compared with hydrody- 
namic simulations for ry/s=0-0.16. 



FIG. 4: ALICE measurements for integrated elliptic flow 
are compared with hydrodynamic simulations for viscosity to 
entropy ratio ?7/s=0-0.16. 



around the critical temperature T = T c [11]. Constant 
■q/s should be understood as space-time averaged r]/s. 
The remaining parameter, central energy density ei at 
the initial time Ti— 0.6 fm is fixed to reproduce ALICE 
measurements for charged particles multiplicity in 0-5% 
Pb+Pb collisions, (^)e* = 1601 ±60 [13]. We sim- 
ulate Pb+Pb collisions and compute the negative pion 
multiplicity. Resonance contributions are included. As- 
suming that pions constitue «85% of the charged parti- 
cles, negative pion multiplicity is multiplied by the factor 
2 x 1.15 to compare with experimental charged particles 
multiplicity. Irrespective of fluid viscosity, in the hydro- 
dynamic model, experimental multiplicity can be fitted 
by changing the initial energy density, more viscous fluid 
requiring less energy density. Multiplicity is proportional 
to final state entropy. In viscous fluid, during evolution, 
entropy is generated. Entropy generation increases with 
viscosity and more viscous fluid require less initial en- 
ergy density to reach a fixed final state entropy. In ta- 
blelH central energy density and temperature of the fluid 
required to reproduce charged particles multiplicity in 0- 
5% collisions are listed. ALICE collaboration measured 
multiplicity quite accurately and initial energy density of 
QGP fluid in ^^^=2.76 Pb+Pb collisions can be esti- 
mated accurately, within 5%. In tableU the bracketed 
values are estimated central energy density and temper- 
ature in -^5^^=200 GeV Au+Au collisions @. From 
RHIC to LHC, collision energy increases by a factor of 
~14 but central energy density of the fluid increase by a 
much lower factor ^4-5. 



V. RESULTS 



participant pairs are compared with the ALICE measure- 
ments. Even though, fluid was initialised to reproduce 
experimental multiplicity in 0-5% collisions, in periph- 
eral collisions, viscous fluid produces more particles than 
an ideal fluid. The reason is understood. Viscous effects 
depend on the velocity gradients. Velocity gradients are 
more in a peripheral collisions than in a central colli- 
sions, so is the viscous effect. To be quantitative about 
the fit to the data in ideal and viscous hydrodynamics, 
we have computed x 2 values for the fits, they are listed 
in tableJIIJ Best fit to the data is obtained in the ideal 
fluid approximation. In viscous evolution \ 2 values in- 
crease by a factor of 2-3.5. We can conclude that ALICE 
charged particles multiplicity data in y/s NN =2.76 TeV 
Pb+Pb collisions do not demand any viscosity. 



TABLE II: X 2 /N values for the fit to the ALICE data for 
different values of viscosity to entropy ratio. 





X 2 /N 


Data set 


T)/S = 0.0 


77/s=0.08 


77/s=0.12 


»7/s=0.16 


charged particles 
multiplicity 


1.46 


3.12 


4.44 


5.24 


0-5% 
pr spectra 


4.38 


7.14 


9.78 


14.84 


V2 integrated 


9.68 


3.02 


12.46 


41.44 


v 2 : 10-20% 


7.07 


4.69 


7.68 


15.14 


v 2 : 20-30% 


23.64 


8.32 


8.46 


18.97 


v 2 : 30-40% 


44.60 


13.14 


9.89 


24.47 


v 2 : 40-50% 


33.62 


9.86 


4.18 


9.00 



A. Centrality dependence of charged particles 
multiplicity 

With all the model parameters fixed we can study the 
centrality dependence of simulated charged particles mul- 
tiplicity in ^/s NN =2.76 TeV Pb+Pb collisions. In FigE] 
simulation results for charged particles multiplicity per 



B. Charged particles pr spectra in 0-5% collision 

ALICE measurements for charged particles px spec- 
tra in 0-5% collisions are compared with hydrodynam- 
ical simulations in Fig|3J We have initialised the fluid 
to reproduce charged particles multiplicity in 0-5% col- 
lision. Such a initialisation do reproduces experimental 
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P T (GeV) 1 D. Differential elliptic flow in 10-50% collision 



FIG. 5: In four panels, ALICE measurements for charged 
particles elliptic flow in 10-20%, 20-30%, 30-40% and 40-50% 
Pb+Pb collisions are shown. The solid, short dashed, dashed 
dotted and long dashed lines are hydrodynamic model simu- 
lation for elliptic flow for fluid viscosity rj/s=Q, 0.08, 0.12 and 
0.16 respectively. 



Pt spectra reasonably well. Effect of viscosity on spectra 
is also evident. Compared to ideal fluid, in viscous evo- 
lution, high pt production is enhanced. In tablelU \ 2 
values for the fits are listed. Here also, we can conclude 
that ALICE data on charged particles pr spectra in 0-5% 
collisions, do not demand any viscosity. 



C. Centrality dependence of integrated elliptic flow 

ALICE collaboration employed various methods, e.g. 
2 and 4 particle cumulant, q-distribution, Lee- Yang zero 
method etc, to measure charged particles elliptic flow 
in Pb+Pb collisions [l2j]. Flow in 4-particle cumulant 
method, Lee- Yang zero method and q-distribution is con- 
sistent with each other. 2-particle cumulant method how- 
ever results in ~15% higher flow. In Fig |4] ALICE mea- 
surements for integrated elliptic flow in 4-particle cumu- 
lant method are compared with hydrodynamical simula- 
tions for elliptic flow. Ellitic flow depend sensitively on 
viscosity, flow reducing with increasing viscosity. Reduc- 
tiuon is more in peripheral than in central collisions. As 
argued previously, viscous effects are more in peripheral 
than in central collisions. For high viscosity, ?7/s=0.16, 
in very peripheral collisions, experimental centrality de- 
pendence is not reproduced. Instead of decreasing V2 
increases, x 2 /N values for the fits to the data are noted 
in table |nl Unlike the centrality dependence of charged 
particles multiplicity or px spectra, best fit to the ALICE 
data on integrated elliptic flow is obtained for r;/,s=0.08. 



In Fig 03 ALICE measurements of elliptic flow «2{4}, 
in 4-particle cumulant method, in 10-20% ,20-30%, 30- 
40% and 40-50% collision centralities are shown. ALICE 
collaboration measured elliptic flow upto pt ~5 GeV. 
Differential w 2 reaches a maximum «0.2 at Pt=3 GeV. 
In FigfSl we have shown the measurements only upto 
Pt=3 GeV. Hydrodynamical model are not well suited 
for large pr- Sources other than thermal also contribute 
to the particle production at large pr- The continuous, 
small dashed, dashed dot, and medium dashed lines in 
FigJS] are simulated elliptic flow for pions in hydrody- 
namic evolution of fluid with viscosity to entropy density 
ratio ?y/s=0.0,0.08,0.12 and 0.16 respectively. 

Unlike the centrality dependence of multiplicity or the 
Pt spectra in 0-5% collisions, differential elliptic flow is 
not best explained in ideal fluid evolution. x 2 /N val- 
ues for the fits, in each collision centralities, are noted 
in tableHU Experimental data only in the pt range, 
0.25 < pt < 3.25 GeV are included in the x 2 analy- 
sis. x 2 /N values indicate that the data are better ex- 
plained in viscous evolution than in ideal fluid evolution. 
It is also apparent that more peripheral collisions demand 
more viscosity. For example, in 10-20% collision, \ 2 is 
minimised for the KSS bound of viscosity to entropy ra- 
tio, 7]/s «0.08. In more peripheral 30-40% collision, data 
demand r\j s ssO.12. The result is interesting. In RHIC 
energy collisions also, it was seen that elliptic flow in pe- 
ripheral collisions demand higher viscosity than in central 
collisions [33]. As discussed in detail in [33|, centrality 
dependence of viscosity can be understood in terms of 
continuously increasing role of hadronic phase of the fluid 
in peripheral collisions. The space-time averaged r\j s has 
contributions, both from the QGP and hadronic phase. 
Compared to a central collisons, in a peripehral collision, 
contribution of hadronic phase to the space-time aver- 
aged rj/s increases, while that of QGP decreases. 

In the present analysis different r\j s are indicated by 
different data sets. For example, ALICE data on mul- 
tiplicity and 0-5% pt spectra are best explained with 
■q/s «0, the integrated vi data with r\j s ssO.08, the dif- 
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ferential V2 data in 10-50% collision with ry/s=0.08-0.12. 
To obtain r\j s which best explain the ALICE data, x 2 /N 
values for the combined data set, (i) centrality depen- 
dence of multiplicity, (ii) centrality dependence of inte- 
grated elliptic flow, (iii) pr spectra in 0-5% collisions and 
(iv) differential elliptic flow in 10-50% collisions, are com- 
puted. They are shown in FigjBJ x 2 /N values exhibit 
a broad minima, indicating that the data are not very 
sensitive to viscosity to entropy ratio. The solid line in 
Fig|5] is a parabolic fit to the x 2 /iV values. From the 
minimum of the fit we extract, rj/s = 0.06 ± 0.02. As it 
was at RHIC enery collisions, at LHC also, experimen- 
tal data are consistent with nearly ideal fluid. We may 
note that the extracted value of rj/s at LHC energy is 
compare well with rjls extracted in other hydrodynamic 
model simulation (34. |35|. However, hybrid models, i.e. 
hydrodynamics coupled to hadron cascade models indi- 
cate higher viscosity [H, H3| • 

In the present analysis we have used some specific 
initial conditions, e.g. initial time t;=0.6 fm, initial 
zero fluid velocity, hard scattering fraction 90%, boost- 
invariant shear stress tensor etc. All possible initial con- 
ditions are not explored. In systematic uncertainty 
in rj/s due to various uncertainties in initial conditions 
was estimated as large as 175%. The presently extracted 
value rj/s = 0.06 ± 0.02 will be even more uncertain, if 



all possible initial conditions are accounted for. 



VI. SUMMARY AND CONCLUSION 

In Israel-Stewart's second order theory of hydrody- 
namics, we have analysed ALICE data on the centrality 
dependence of charged particles multiplicity, integrated 
and differential elliptic flow, charged particles pr spec- 
tra in 0-5% collisions. ALICE data for charged paricles 
multiplicity or 0-5% px spectra do not demand any vis- 
cosity, data are best explained in ideal fluid approxima- 
tion. On the other hand, ALICE data on integrated or 
differential elliptic flow demand small viscosity. From a 
simulataneous analysis of all the data sets we obtain vis- 
cosity to entropy ratio r]/s = 0.06 ±0.02, in y/s NN —2.76 
TeV Pb+Pb collisions. The extracted values of viscosity 
is rather similar to the value extracted at RHIC energy, 
even though the fluid at LHC is at higher temperature 
than at RHIC. It appear that both at RHIC and LHC, 
nearly perfect fluid is produced. 

Acknowledgement: Authors would like to thank Dr. 
Bedangadas Mohanty for providing the ALICE data on 
charged particle elliptic flow for different centralities. 
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